Fronts, Domain Growth and Dynamical Scaling in a d = 1 non-Potential System 
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We present a study of dynamical scaling and front motion in a one dimensional system that 
describes Rayleigh-Benard convection in a rotating cell. We use a model of three competing modes 
proposed by Busse and Heikes to which spatial dependent terms have been added. As long as the 
angular velocity is different from zero, there is no known Lyapunov potential for the dynamics of 
the system. As a consequence the system follows a non-relaxational dynamics and the asymptotic 
state can not be associated with a final equilibrium state. When the rotation angular velocity is 
greater than some critical value, the system undergoes the Kiippers-Lortz instability leading to a 
time dependent chaotic dynamics and there is no coarsening beyond this instability. We have focused 
on the transient dynamics below this instability, where the dynamics is still non-relaxational. In 
this regime the dynamics is governed by a non-relaxational motion of fronts separating dynamically 
equivalent homogeneous states. We classify the families of fronts that occur in the dynamics, and 
calculate their shape and velocity. We have found that a scaling description of the coarsening process 
is still valid as in the potential case. The growth law is nearly logarithmic with time for short times 
1 and becomes linear after a crossover, whose width is determined by the strength of the non-potential 

terms. 
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I. INTRODUCTION 

' A noteworthy result in the study of non-equilibrium statistical mechanics is the existence of dynamical scaling 
during the coarsening process in which a system approaches equilibrium after undergoing a phase transition 
Dynamical scaling reflects that domain growth is self-similar with a single time dependent characteristic length. In 
the simplest case of a relaxational dynamics for a scalar order parameter, which models, for example, an order-disorder 
I ■ transition (model A in the taxonomy of S), domains of two equivalent phases grow locally from an unstable state, 
and the approach to a final equilibrium state is dominated by interface motion. For spatial dimension d > 1, the 
mechanism for domain growth is curvature driven interface motion. This leads to a characteristic length growing as 
R ~ t 1 / 2 . This type of phenomena has been studied in a large variety of systems which share the common feature that 
the final state of the dynamics is a state of thermodynamic equilibrium which minimizes a free energy or "potential" 
of the problem. Transient dynamics might also include additional processes beyond pure relaxation in that potential 
k> , J|,[|], but a measure of relative stability between stationary states is guaranteed by the existence of the potential. A 
more genuine non-equilibrium dynamics occurs when such a potential does not exist. A natural question, which we 

[ address in this paper, is the existence of dynamical scaling in the approach to a final stationary state which does not 
follow the minimization of a potential, while this transient dynamics from an unstable state involves the formation 
of spatial domains. The question of dynamical scaling can also be addressed for Hamiltonian dynamics ||, which is 
the extreme opposite situation to that of dissipative relaxational dynamics in a potential. A general non-equilibrium 
situation will, in general, have contributions from both types of dynamics [Q,|5| . 

The motion of an interface between two linearly stable solutions of a dynamical system was long ago proposed 
as a measure of relative stability for a non-potential system and the motion of interfaces, domain walls or front 
solutions has been studied in a number of non-potential systems ]^-pd| . It is known that a domain wall between two 
equivalent states with different broken symmetry can move in d = 1 in either direction due to non-potential dynamics 
Likewise, non-potential dynamics can stabilize front solutions which, in the potential limit, would move from a 
globally stable into a metastable state It is the purpose of this paper to study the consequences of interface 

motion driven by non-potential dynamics on coarsening processes and dynamical scaling, which have been by and 
large not considered. 
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A physical prototype situation in which the question posed above can be studied is Rayleigh-Benard convection in 
a rotating cell |l2||, w hich in many aspects is mathematically equivalent to a biological model of competition between 
three species Jl3| , |l4| . Beyond a threshold value for the parameter, 8, measuring the strength of non-potential terms 
(related to the rotation speed in the case of Rayleigh-Benard convection) an instability to a time dependent dynamics 
occurs (Kiippers-Lortz instability |l5|,[l6|] ) . Below this instability, but still taking the system beyond the Rayleigh- 
Benard convective instability, locally ordered domains, associated with different orientations of the convective rolls 
emerge. The subsequent coarsening process seems to be stopped by the Kiippers-Lortz instability in d = 2 |17|^8| . 
However, below the instability to a time dependent state, three preferred orientations exist and the motion of interfaces 
separating them is subject to non-potential dynamics which will affect domain growth. 

With this motivation in mind, and as a first step towards the understanding of the problem of domain growth and 
dynamical scaling in this type of systems, we have considered a d = 1 model for three competing real non-conserved 
order parameters with the non-potential dynamics of [ fulfill and spatial inhomogeneities determined by short range 
self- interactions. We have chosen this system because dynamical scaling in d = 1 potential systems is a special case 
for which well established results are available for the simplest case of a scalar non-conserved order parameter 
with short range interactions, a scaling solution is known with a logarithmic growth law for the typical domain size, 
R r~j logi [|l9|,g0). This regime follows an early time regime of domain formation with a growth law R ~ t 1 / 2 pi) . 
The logarithmic domain growth has its origin in the interactions between domain walls ||22| . The velocity of domain 
wall motion in these circumstances can be calculated by a perturbation analysis |23| . We study how these results are 
modified by non-potential dynamics. We find that dynamical scaling still holds, but with a crossover between two well 
defined regimes characterized by a logarithmic and linear domain growth law respectively. The two growth laws can 
be traced back to the two mechanisms that determine domain wall motion. The first one is the interaction between 
domain walls as in the potential case. The second one is due to the fact that the non-potential dynamics causes 
that isolated individual fronts move with finite velocity. In a multi-front configuration this provides and additional 
coarsening mechanism in which fronts moving in opposite directions annihilate each other. The crossover time between 
the two dominant mechanisms described depends on the strength of non-potential terms. When these become large 
enough, the logarithmic regime is pushed to just the very early times. In this case finite size effects become also 
important since very large domains emerge rather fast. 

The outline of the paper is as follows: in section || we introduce the non-potential model and describe its homoge- 
neous stationary solutions. In section III we clasify non- homogeneous front solutions and compute their velocity. We 
also analyze the interactions between two fronts including the non-potential effects. In section IV we discuss the issue 
of dynamical scaling and present numerical simulations that show the validity of the dynamical scaling description 
when non-potential terms are present. Finally, in section [v] we end with some conclusions and an outlook. 



II. THEORETICAL MODEL 



We base our theoretical approach in a three mode model first proposed in the context of fluid dynamics by Busse 
and Heikes Jig ). In this model, the (real) amplitudes of the three selected modes corresponding to three different 
orientations of the convection rolls, A\, A2, A3, follow the evolution equations: 
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(2.1) 

A similar set of equations was proposed earlier to model the behavior of three competing biological species [ fl3| . In this 
case, Ai, A 2 and A3 stand for the population number of each species. For the fluid case, 8 is related to the rotation 
speed such that <5 — is the non-rotating case; the parameter rj is related to other fluid parameters. The analysis of 
|l6f and |l3f| shows that, for a certain range of the parameters r\ and 8 (see next section) there are no homogeneous 
stable solutions and the dynamics tends asymptotically to a sequence of alternations of the three modes as shown in 
experiments. An unwanted feature of the previous model is that the alternation time is not constant, but increases 
with time, contrary to experiments where an approximately constant period is observed. Although Busse and Heikes 
proposed that the addition of small noise could stabilize the period P]24|] , an alternative explanation considered the 
addition of spatially dependent terms to the previous equations. Although the symmetries that must satisfy the 
amplitude equations would imply that the spatially dependent terms should be of an specific form p5|-p7| , it has been 
shown in Jl7] that these can be further simplified without altering the essentials of the problem. We choose in this 
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work the simplest diffusive form for the spatial depending terms, namely: 



dtAi = d 2 x A x + A 1 (l-A 2 -(r ) + 8) A\ - (r, - 8) A 2 ) 

d t A 2 = dlA 2 + A 2 (1-A 2 -( V + 8) A 2 -{ v - 8) A\) (2.2) 
d t A 3 = d 2 x A 3 + As (1 - A 2 - (r, + 8) A\ - (r, - 8) A\) 

which form the basis of our subsequent analysis. A similar set of equations for the modulus square of the amplitudes 
is introduced in reference fl4|| for some particular values of the parameters T), 8. Notice that the system is invariant 
under the following transformations: 

a) x —>■ x + xq, t —> t + to (spatio-temporal translation symmetry) 

b) A\ — ► A 2 , A 2 — * A3, A 3 — ► A\ (cyclic permutation symmetry) 

c) Ai Aj, 8 — ► — <5, where Aj, Aj are any two different amplitudes. 

The previous analysis shows that A%, A 2 and A3 are "equivalent" (from the dynamical point of view) variables. 
It is possible to split the dynamical equations into a potential and a non-potential contributions: 

d t A i = -—-6-f i , i = 1,2,3 (2.3) 
oAi 

where the potential function J- is given by: 



F[A l7 A 2 ,A 3 } = - f dx(- [(d x A x ) 2 + (d x A 2 ) 2 + (d x A 3 ) 2 ] + (Aj+Al + Aj) 

J , 1 J , (2.4) 

--{A\ +A* + Ai)-± V (A 2 A 2 + A 2 A 2 + A 2 A 2 ) 



and the non-potential terms are: 
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In the case 6 = the dynamical flow is of type relaxational gradient |4|,|5j], that is, there exists a Lyapunov functional 
(JF) that monotonically decreases in time. When 8 ^ 0, however, one cannot find generally such a functional and we 
say that the system is non-potential. 

There exist two kinds of homogeneous solutions which are stable in some region of the parameter space spanned 
by 77 and 8. These are three "roll" solutions Aj = 1, Aj = (j ^ i), i = 1,2,3, and one "hexagon" solution 
A\ = A 2 = A3 = l/v 7 ! + 2?7 (this solution requires 77 > —1/2). A linear stability diagram of these solutions is shown 
in figure [l]. We have focused on the region labeled with the letter 'R' for rotation speeds below the Kuppers-Lortz 
instability (\S\ < rj — 1) and where the rolls are the stable solutions. In this region we expect the formation of domain 
walls connecting homogeneous stable roll solutions. 

III. FRONT SOLUTIONS 
A. Isolated fronts 

In the context of the present study, fronts or domain walls are defects that connect two stable homogeneous 
solutions. Fronts in one dimension are usually termed kinks and we will often refer to them in this way. We focus 



on the spatial-dependent stationary solutions of (2.2) with 8 = 0. The stable kink solutions are such that one of the 
three amplitudes, say A^, satisfying the boundary conditions Ak(x — » ±00) = is zero everywhere. In order to study 
the dynamics of the non-potential kinks, we first consider the kinks associated with the stationary potential problem 
and then we treat the non-potential terms as a perturbation. The two non-vanishing stationary amplitudes Aj and 
Aj are, for 8 = 0, solutions of: 

d 2 x A z = -A + Af + V AA 2 

• a; ■ 
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with boundary con ditions Aj(— oo) = A,(+oo) — and Ai(+oo) = Aj(— oo) = 1. 

The system (3T) may be considered to represent the two dimensional motion of a Newtonian particle of unit mass 
(a; -> t, A l -> X, -> Y") under the action of a force with potential function V(X, Y) = \(X 2 + Y 2 ) - \(X i + 
Y" 4 ) — \r\X 2 Y 2 . This function has two maxima in mo = {Ai = l,Aj = 0} and mi = {Ai = l,Aj = 1}. It is clear 
that there exists a unique trajectory (allowed by the dynamics) along which a particle located in mo (mi) can reach 
mi (mo). The kink profile corresponds to the variation in time of the particle coordinates (X(t),Y(t)) when it moves 
between the two maxima p8|| . 

An explicit analytical solution can be found in two particular cases p9[. First, when < T) — 1 <C 1, we have: 



A°{x) = r(x) [l + exp (2^-1(1 - a? )) 



-1/2 



A°(a?) = r(x)e x 1 + exp ( 2y/rj - l(x - x ) 



-1/2 



(3.2) 



with r(x) = 1 + [r) — l)R(x), R(x) — 0(1). Secondly, when i] 
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3 it is possible to obtain exact analytical solutions: 

(3.3) 
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In both cases Xo is arbitrary but fixed. From these solutions it is clear that the spatial scale over which A? and A° 
vary is of order l/y/rj — 1. 



The three roll solutions are equivalent and they yield the same value for the Lyapunov functional (2.4). Therefore, 
we expect kinks not to move in the potential problem (S = 0). We now ask about the persistence of these kink 
solutions when 6 is different from zero. For this we will use singular perturbation theory. Let us assume S to be small, 
say of order e, and look for a solution of (p.2[) (with A^{x) = 0) of the form: 



(3.4) 



Ai(x) = A a t (x - s(t)) +eA\{x - s{t)) + 0{e 2 ) 
Ajix) = A°Jx - s(t)) +eA){x - s{t)) + 0(e 2 ) 



where A®(x) and Aj(x) are solutions of (3.1). Substituting into (2J2) and matching the terms of the same order in e, 
we find, at order 0(e°): 
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and at order 0(e 1 ): 
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The solvability condition for the existence of a solution (A](x), Aj(x)) reads 

(* f , a')=0 



(3.7) 



where (•, •) is a scalar product in i 2 (M) defined by (/, g) = j x ' oo dxf(x)*g(x) and & belongs to the null space of the 
auto-adjoint operator C. Because of the translational invariance, C has a zero eigenvalue so that its kernel is not 
empty. The associated eigenvector is generally known to be 
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This is immediately seen taking, for example, the derivative of (3.5) with respect to x. Eq. (3.7) can now be explicitly 
evaluated. From this equation , the solitary kink velocity in the non-potential case is obtained at leading order: 



r dxASA AA <d x AS-J8d 9 A<l) 

v ^ ^ dtS = S r° a Wh \ l^lm 40,21 + °( S ) ( 3J ) 
J_ 00 dx [{d x A(Y + (d x A°f\ 

Therefore, in the non-potential case, the kink moves despite connecting states associated with the same value of 
the Lyapunov potential of the equilibrium problem, as already known for other problems . For the particular case 
of t] = 3 for which an analytical result is available for the kink profile, (eq. (3.3)), an explicit result is obtained for 
the solitary kink velocity, namely v(S) = 6^/2/4 



The expression (3.9) gives not only the magnitude of the velocity but also the direction of motion, which is related 
to the sign of v. First, we note that the velocity is at leading order proportional to 5, so the direction of the motion 
depends upon the sign of S. To illustrate how ( |3.9D determines the direction, let us consider for example a kink with 
boundary conditions: oo) = A, (+oo) = a nd Aj(+oo) = Aj(— oo) = 1; A® (x) and Aj(x) are such that d x A® > 
and d x Aj < 0. In this case the numerator of ( |3.9p is positive and v has the sign of S. A positive (negative) value 
of v corresponds to a kink moving to the right (left). In figure || we show a classification of the six possible types of 
isolated kinks and their direction of motion. Three of them move in one direction and the other three in opposite 
direction. 

We have checked numerically the domain of validity of the perturbative result fl3.9| ) (see figure ||) . To check fl3.9| ) 
we either use the analytical result of the kink profile A® for r\ = 3, or, more generally, the k ink profile A® obtained 



numerically. For a value of r\ = 3.5, we see that the perturbative result to first order in 8 (3.S) turns out to be in 
good agreement with the numerical results approximately for values of 8 < 1.5. Of course this upper limit of validity 
depends on 77 in such way that it gets bigger as rj is larger. Above this limit the linear relation between v and 8 is no 
longer valid and one needs to compute further corrections in terms of successive powers of 8. 



B. Multifront configurations 



To study transient dynamics and domain growth we consider random initial conditions of small amplitude around 
the unstable solution A\ = A2 = A3 = 0. In this situation a multifront pattern emerges rather than a solitary kink. 
In order to study dynamical scaling, we are interested in the late stage of this dynamics, once well-defined domains 
have been formed. 

In a potential system governed by a non-conserved scalar order parameter with short range interactions, as it 
the case with 5 — 0, late time dynamics can be explained in terms of the interaction (and further annihilation) 
among adjacent kinks [pT |l9[j 20|j2^| . An isolated kink is stable. The interacting force between kinks turns out to be 
proportional to exp(— ad) ||22|| , where a is some positive constant related to the interface width and other system 
parameters, and d is the distance between two adjacent kinks. This interaction among kinks leads to a growth law 
for the characteristic length that depends on time logarithmically ||l9| , p0| . The force is attractive and leads to kink 
annihilation. The process occurs in such way that the domain that vanishes first is the smallest one. Kink annihilation 
occurs in a very small time scale. In fact, the hypothesis of "instantaneous annihilation" has been found to be a good 
assumption [[L9| . Kink annihilation induces domain coarsening leading to a final state with an homogeneous roll 
solution filling up the whole system. 

When 8 is different from zero the long stage dynamics should still be explained in terms of moving fronts which 
annihilate each other. But now two very distinct competing physical phenomena come into play. On the one hand, 
there is the aforementioned kink interaction. On the other hand, we have the kink motion driven by non-potential 
effects. In this case, we do not expect the growth law to be logarithmic, at least in the regime where the non- 
potential effects (the strength of which is measured by 8) are important. In figure (jj) we show some snapshots 
corresponding to a typical run of the temporal evolution of the system (we use periodic boundary conditions). The 
first snapshot corresponds to an early stage during which domains are forming. Once formed, kinks move in such a 
way that annihilations of contrapropagating adjacent kinks leads to coarsening. Eventually, as corresponding to the 
last snapshot, the system may be in a state with a group of kinks moving all in the same direction. These will interact 
among them (with a interaction force that varies logarithmically with the interkink distance) until extinction. 

We have performed a perturbation analysis of domain growth in the simplest example of a single domain bounded 
by two moving domain walls. A differential equation for the domain size s(t) can be obtained in the general case (see 
appendix). For 77 = 3 it adopts the simple form: 

d t s(t) = 2v(S) - 24\/2e- v5s W (3.10) 
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where v(S) is the solitary kink velocity. This expression is obtained in the "dilute-defect gas approximation", that 
is, when the width of the fronts is much smaller than the distance between them. The first term in the right hand 
of ( |3.10| ) can be either negative or positive and represents the contribution to the variation of the domain size owing to 
non-potential effects. The second term is related to the interacting force between the kinks and it is always negative 
(attractive force) so that it tends to shrink the domain. If both terms are negative, the kinks will annihilate each 
other. Otherwise, when the first term is positive, the two effects act in opposite directions. In fact, given an initial 
size of the domain so, it is possible to find a value 5 = 5 C for which the domain neither shrinks nor grows; in this 
case, the initial domain would not evolve in time being a stationary solution. For values S > 5 C the domain would 
get wider whereas for 6 < 6 C it would shrink. Note that in the 5 = case (potential regime), the isolated domain 
always collapses, but this can be stopped with a suitable strength of the non-potential terms. Furthermore, given a 
fixed value of 5, if Sq is large enough, the dominant term responsible of the kink motion is the one associated with 
v(5). In this case the fronts move at a constant velocity leading to a variation of the domain size linear with time. 
On the other hand, if sq is small enough, kink interaction will be the dominant effect and the single domain size will 
collapse logarithmically with time. This picture of the size dynamics of a single domain also explains basically what 
happens when more domains (and a non-vanishing third amplitude) coexist. It gives a useful understanding of the 
characteristic growth laws obtained from a statistical analysis in the next section. 



IV. DOMAIN GROWTH AND SCALING 



In this section we focus on the scaling properties of the system (2.2) in a late stage of the dynamics, namely when 
well-defined domains have formed. The scaling hypothesis states that there exists a single characteristic length scale 
R(t) such that the domain structure is, in a statistical sense, independent of time when lengths are scaled by R(t). 
We will refer to the time dependence of the scale length as the growth law of the system. It has been found that the 
scaling hypothesis holds in a great variety of potential systems. The system under study here gives us the opportunity 
of answering the question of whether a non-potential dynamics satisfies dynamical scaling. 

Two magnitudes frequently used to study domain growth and scaling properties for a scalar field ^(x, t) (for 
instance, one of the three amplitudes in equations (2.2)) are the equal time correlation function 

C(r, t) = *( x + r, *)y (4- 1 ) 

and its Fourier transform, the equal time structure factor 

S(k,t) = ^*(M)*(-M)^ (4.2) 

where the angular brackets indicate an average over initial conditions ("runs"). If a single characteristic length exists, 
according to the scaling hypothesis, the pair correlation function and the structure factor must have the following 
scaling forms in a d-dimensional system: 

C(r,t) = f(r/R(t)) (4.3) 
S(k,t) = R(t) d f(kR(t)) (4.4) 
The function / is called the scaling fun ction. To check numerically the validity of the previous scaling laws, we have 



integrated the system of equations (2.2) using a finite difference method for both, the spatial and temporal derivatives. 
In the simulations we have taken a constant value for 77, namely 77 = 3.5, and we have varied S from 6 = (potential 
case) to 5 = 0.1 (a value below the KL instability threshold). We have used periodic boundary conditions and have 
averaged our results over 100-500 runs and used system sizes ranging from L = 250 to L = 1000. To study domain 
structure, we consider the correlation function of one of the three amplitudes. We use the correlation function better 
than the structure factor because of the large fluctuations of the structure factor at small wave numbers. A typical 
length scale associated with the average domain size can be defined in several ways. Specifically, we have determined it 
by computing the value of r for which C(r, t) is half its value at the origin at time t, that is C(R(t),t) = |C(0, t). The 
calculation has been performed by fitting the four points of C(r, t) closest to C(0, t) to a cubic polynomial. Another 
typical length, Ri (t) can be evaluated directly as the system size divided by the number of kinks. We have verified 
that the quotient R±(t) / R(t) remains nearly constant, as expected, when a single characteristic length dominates the 
problem. 
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A. Growth Law 



We consider first the potential case 5 = 0: In figured we show that the domain size follows the expected logarithmic 
behavior. The attractive interaction among the kinks leads to a very long transient before the system reaches its final 
state which corresponds to one roll solution filling up the whole system. 

In the non-potential case the domain size R(t) is shown in figure || for a system of size 500 and 5 = 10~ 3 . For the 



earliest times, when the kinks are very close to each other, and according to the discussion in section III B we expect 
the interaction terms to be the dominant ones (as long as S is small enough). This leads to a logarithmic growth of 
R(t) as observed in region (a) of figure ^. Due to coarsening the characteristic domain size becomes larger and the 
domain wall interaction becomes weaker as the time increases. For longer times the non-potential effects dominate 
with respect to wall interaction. In this regime we can consider each domain wall to move at a constant velocity This 
gives rise to a linear behavior of R(t) with time (region (7)). Between regions (a) and (7) there exists a crossover 
(region (/?)) for which the weights of both effects (interaction and non-potential) in driving the domain wall motion 
are of the same order. Finally, at very late times finite size effects come into play (region (6)): the domain size 
saturates to a constant value and the number of domain walls is too small to make good statistics. 

In the regime for which the interaction effects are the dominant ones, we can give a simple explanation of the linear 
growth law observed for R(t) . Statistically speaking, there will be the same number of kinks moving to the right and 
to the left. As a matter of fact, in an appropriate reference frame, the system can be seen as composed of motionless 
kinks (type I) and kinks moving at a velocity of 2v in one fixed direction (type 2) . If we call Ni(t) and N 2 (t) the 
average number of kinks of both types at time t, the number of kinks of, say type 1, at time t + dt will be: 

Ni(t + dt) = Ni(t) - Ni(t) v dt x (4.5) 

±j 

where the second term in the right hand side represents the number of kinks disappeared in dt by annihilation and L 
stands for the system size. The important point is that N\(t) = N^it) = N(t) (remember we are dealing with averaged 



quantities), so that (4.5) transforms into 



« = 1 -| W(t) * (4.6) 

The integration of the previous equation gives N(t) = [(v/L)t + N^ 1 }^ 1 ~ t -1 , so that the average inter-kink distance 
R(t) ~ iV(i) -1 ~ t is linear with time. 

When 5 is large enough the initial kink annihilation is so fast that the regions (a) and (f3) in the R(t) plot can 
hardly be observed in the numerical integration. In this case of large non-potential effects, a linear growth law is 
observed for the shortest times and it continues until finite size effects show up clearly (see figure |?]). Finite size effects 
are seen to occur for t > 200 for S = 0.1. The bigger 5, the sooner finite size effects appear. 

At long times the system will consist of an homogeneous roll state or a group of kinks moving either to the right 
or to the left pl[ | . Note that the periodic boundary conditions impose constraints about the number of such moving 
kinks. To be precise, the number of kinks moving in a fixed direction must be multiple of three. We can form a 
subgroup of three kinks moving in the same direction by joining those appearing in each row of figure The 
moving kinks will carry on interacting among them (logarithmically) until eventually they all will disappear. In this 
situation we expect the growth law to be logarithmic with time but one of such groups is composed typically of three, 
six or rarely nine kinks, a number too small as to make good statistics. 



B. Scaling function 



We now address the question of the validity of the dynamical scaling hypothesis (4.3). For this purpose, we have 
plotted the equal time correlation function C(r, i) versus the scaled length r/R(t) for several times. Figure || shows 
the scaling function in the potential case. In figures [)]and [h] we show the correlation functions for several times before 
and after scaling the system length. Our results show that the correlation functions follow a single profile when the 
length is scaled with the characteristic domain sizes obtained above. We therefore conclude that a scaling description 
of the system is also valid as in the potential case, but now with a non-potential dynamics. The upper limit of the 
time interval during which there is scaling is determined by the appearance of finite size effects. The range of values of 
S for which there is scaling in a quite large time interval is rather small. For values of 6 even of a few tenths, the finite 
size effects show up for very short times. Moreover the fluctuations in the scaling function grow as S increases. For 
these reasons, we have not been able to obtain a conclusive comparison among scaling functions for different values 
of 5, although their shapes appear to be rather insensitive to the value of S. 
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V. CONCLUSIONS 



We have studied domain growth and dynamical scaling in a non-potential coarsening process in one dimension. The 
model considered is motivated by the phenomenon of Rayleigh-Benard convection in a rotating cell. We have focused 
on the region below the Kuppers-Lortz instability point, where the dynamics is still non-potential and the system 
shows coarsening. A solitary kink moves at a constant velocity due to the non-potential dynamics. When there are 
several kinks present in the system, these move due to both, domain wall interaction and non-potential effects. In any 
case the dynamics is governed by motion of interfaces. This motion is such that kinks moving in opposite directions 
annihilate each other. As a consequence of kink annihilation the average domain size grows in time and the system 
coarsens. When 6 = (potential case) we have shown that, in accordance with general results, the growth law is 
logarithmic with time and that a scaling description of the system dynamics is possible. When S becomes different 
from zero we have found that the scaling hypothesis still holds, as in the potential case, but with a different growth 
law that reflects the non-potential dynamics of the system. For the shortest times, the kink interaction (the only 
effect present in the potential case) is the dominant effect and gives rise to a logarithmic growth law with time. For 
longer times the average inter-kink distance is large enough to make the interaction effects negligible in driving kink 
motion. Therefore each kink moves nearly independently of the others as if it were isolated. In this situation each 
kink moves with a constant velocity leading to a linear growth law. For larger values of 6 the logarithmic region is not 
observed because of the fast annihilation of the domains during the very early times. The two dimensional version of 
this problem is currently under study |32] ] and it exhibits rather different dynamical behavior grossly dominated by 
vertices where three domain walls meet and which have no parallel in one dimensional systems. 



APPENDIX 



We consider an isolated domain bounded by two domain walls associated with amplitudes A\ and A 2 , while A3 = 0. 
When the domain size is much greater than the interface width ("dilute-defect gas approximation"), a reasonable 
ansatz for this solution is 



Ax{x,t) = a(x - r(t)) + b(x -d + r(t)) + w 1 {x,t) 
A 2 (x, t) = b(x - r(t)) +a(x-d + r(t)) - 1 + w 2 (x, t) 



(A.I) 



where r(t) measures the displacement of the kinks, d is the initial domain size (so that the domain size at time t is 
d — 2r(t)), dtr and Wi (i — 1, 2) are assumed to be small corrections of order 5 and dtWi to be negligible with respect 
to Wi. To simplify notation, we use: / = f(x — r(t)), f d = f(x — d + r(t)). The moving fronts a and b satisfy the 
boundary conditions a(oo) = b(— 00) = 0, a(— 00) = 0(00) — 1 and they are solutions of the system (2.2) (with one of 
the amplitudes equal to zero) so that the following equations hold: 



M+{a,b) = M+(b d ,a d ) = M-(b,a) = M-{a d ,b d ) = 
where the action of the operators A4±(-, •) is given by: 

M ± (f,g) =dl± v(S)d x +f-f-( v ± S)fg 2 



(A.2) 



(A.3) 



The parameter v(5) is the fron t velo city as given by eq. ( |3.9| ). 

Introducing the ansatz (Al) into ( |2.2| ) we obtain, at leading order, a linear system of equations for wi and W2'- 



C = 



4> = 



£<j> = <p' 

- 3(a + b d f - n(b + a d - l) 2 -2r)(a + b d )(a d + 6-1) 

-277(0 + b d ) (a d + b - 1) &l + 1 - 3(a d + b - l) 2 - //(a + b d ) 2 

■ 1 _ ( (d x a + d x b d )v(S) + (d x b d - d x a)d t r + KiS + K 2 r] + K3 
Xd x b + d x a d )v{8) + (d x a d - d x b)d t r + K[6 + K' 2 r) + K 3 



(A.4) 



, = 



8 



where the functions Ki(x,t) and K'^x^t) (i = 1,2,3) are given by: 



■Mi 



K x = a(a d - l) 2 + 2b(a + b d )(a d - 1) + b d (l - 2a d + b 2 ) 
K 3 



K' 2 



a(a d - 1) + 2(a d - l)(ab + a d b d + bb d ) 
-2 d 2 x b d + b d (3a 2 + 3ab d + 2b 2 - 2) 
~a(a + 2b d )(a d - 1) + b d {b d - 2ab - bb d ) 
a(a + 2b d )(a d - 1) + b d (2ab + bb d + 2a d b d - b d ) 
-2 d 2 x a d + 3b 2 (a d - 1) + 3b(a d - l) 2 + 2a 3 d - 3a 2 . 



a d 



The solvability condition for the existence of a solution (wi(x,t),w 2 (x,t)) for (A. 4) reads 

(* f , 4>') = o 



(A.5) 



where belongs to the kernel of the auto-adjoint linear differential operator C. We will show below that VP' is 
approximately given by (d x a, d x b) T (here T denotes the transposed vector), where a = a(x — r(t)) and b = b(x — r(t)) 
are the domain wall profiles around x = r(t). 

The first component of the vector is given by: 



(C&)i =£ n d x a + Ci 2 d x b 



d x a + d x a- 3(a + b d ) 2 d x a- r)(a d + b- l) 2 d x a- 2rj(a + b d )(a d + b- l)d x b 



(A.6) 



As long as that the width of the interfaces is much smaller than the domain size (for all times t) 7 we can make the 
following approximations: ab d sa 0, aa d w a, bb d w b d . Moreover, this assumption implies that the product of the 
derivative with respect to x of an amplitude solution centered on x — xq multiplied by another amplitude shifted a 
length of order of the domain size, will be a function which will take values different from zero only in a small region 
around x — xq. By using the approximations 

(a + b d ) 2 d x a w a 2 d x a 



(a d + b-l) 2 d x a^b 2 d x a 
(a + bb)(a d + b — l)d x b i=s abd x b 



we find: 



d x [d 2 x c 



i] b a] 



Taking the derivative of (2.2) with respect to x we find that the right hand side of (A. 
amplitude solutions a — a(x — r(t)) and b = b(x — r(t)) are replaced by its form for 5 = 



(A.7) 



(A.8) 

i) is equal to zero when the 
0. Hence, we conclude that 



(£* T )i = O(S). Likewise, we can prove that (C^')2 = 0(8). Therefore, at lowest order in 5, (d x a, d x b) belongs to 
the kernel of the operator C. 

Now we can calculate the evolution of the domain size s(t) = d — 2r(t) through the solvability condition ( |A.5| ). We 
obtain: 



d t s ^ ±2v(8) 



dx(h a d x a + hb d x b) 

jr Qa dx[(d x ar + (d x br] 



(A.9) 



where the coefficients h a and hb depend upon the amplitude solutions a and b and the non-potential parameter 5. 
The first term of the right hand side of (A.9) represents the rate of change of the domain size due to non-potential 
effects which cause the kinks to move at a constant velocity v (S). The s econ d term is related to kink interaction. In 
the case r\ = 3 we can compute explicitly all the coefficients involved in (A.9) taking advantage of the analytical kink 
profiles at lowest order in 5 (equation ( |3.3| )). Making an expansion in powers of e _v/2s ('\ retaining only the leading 
terms, and provided that S is a small parameter, we obtain: 



d t s = ±— - 24V2e 
V2 



-V2s(t) 



(A.10) 



which is eq. (3.10). 
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FIG. 1. Linear stability diagram of the homogeneous solutions of (2.2). Inside the region labeled with R, rolls are stable 
whereas in the H region, the stable solution is the hexagon. The KL region corresponds to the Kuppers-Lortz instability. 
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FIG. 2. Kinds of fronts and their direction of motion for S > 0. The remaining amplitude for each kink is understood to be 
zero across the interface. For S < the picture is the same but with the arrows interchanged. 
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FIG. 4. Snapshots of the temporal evolution of the system. Parameter values: r\ — 3.5, 5 — 0.5, L — 500. 
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FIG. 5. Time evolution of the characteristic domain size for the potential case S = and L = 1000. The straight line is a 
linear regression fit of points obtained numerically. 
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FIG. 6. (a) Time evolution of the characteristic domain size for 8 = 0.001 and L = 1000. The initial logarithmic growth law 
(region (a)) becomes linear (region (7)) after a crossover (region (/3)). Region (S) is related to finite size effects, (b) Zoom of 
region (a) in the left plot with logarithmic time scale. 
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FIG. 7. Time evolution of the characteristic domain size for S = 0.1 and L = 1000. 
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FIG. 8. Scaling function for the potential case 5 = 0. The plot has been made by over plotting C(r, U) vs. r/R{U) for several 
times from t = 200 up to t = 5000. 
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FIG. 9. (a) Equal time correlation function versus the non-scaled length for S = 0.001, L — 1000 and several different times 
from t=150 to t=15000. (b) Equal time correlation function versus the scaled length. The system parameters and the times 
for each curve are the same as in figure (a). 




FIG. 10. (a) Equal time correlation function versus the non-scaled length for S — 0.1, L = 1000 and several different times 
from t=15 to t=150. (b) Equal time correlation function versus the scaled length. The system parameters and the times for 
each curve are the same as in figure (a). 
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